How to use gamma function in integration